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Abstract 

Hadron multiplicity from W boson is calculated in pQCD. The 
agreement of our theoretical predictions with the LEP data says in 
favor of universality of the QCD evolution in hard processes. 



1 Introduction 

Experiments at LEP and SLAC colliders have shown that the multiple pro- 
duction of hadrons in e + e~ annihilation depend on the mass of the primary 
(anti)quarks which launch the process of the QCD evolution. Let us consider 
a heavy quark induced event, 

eV^QQ^I, (1) 

where Q means a heavy (c or b) quark, and an e + e~ event induced by the 
pair of light quarks: 

e + e-^ll^X. (2) 

Here and in what follows I denotes u, d or s-quarks which are assumed to be 
massless. In both cases, X means a system of final hadrons. 
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Let Nqq(W 2 jItiq) and N U {W 2 ) be the average multiplicities of charged 
hadrons in e + e~ events with heavy (JTJ and light primary quarks ([2]), respec- 
tively. W is the invariant energy of colliding leptons, uiq is the mass of 
the (anti) quark Q. It appeared that differences between the light and heavy 
quark induced multiplicities, 

Sqi = N QQ - N a , (3) 

become independent of the collision energy W, but depend only on the heavy 
quark mass tuq. QCD calculations describe the phenomenon quite well p]-[3] 
(see also jl]). 

The QCD calculations of the hadron multiplicities in e + e~ events, in 
particular, hadron multiplicity from Z boson, are in a good agreement with 
the data. Their energy dependence is defined by the QCD evolution of the 
parton showers. The aim of the present paper is to calculate the hadron 
multiplicity from the W boson in pQCD, and thus to check once more the 
universality of the QCD evolution in hard processes. 



2 Hadron multiplicity from the W boson 

From now on we will consider the multiple hadron production in e + e~ anni- 
hilation mediated by production of a pair of W bosons. We will refer to e + e~ 
event as the heavy quark event if one of the W bosons produces a lepton pair, 
while the other decays into hadrons via charm production, for instance: 

e + e-^W + W~ } W + ^cs^X, W~ -»• fi'u„, . (4) 

In the light quark event final hadrons are fragments of the light quark- 
antiquark pair produced by one of the W bosons: 

e + e~ -> W + W~, W + -+IV -> X, W~ -> W~ -> y~v^ . (5) 

Let us define the average multiplicities of charged hadrons in the above 
mentioned processes as Ngi(W 2 ,m 2 ) and Ni(W 2 ), where m is the mass of 
the heavy (charm) quark[j Our main goal is to establish a relation between 
these two multiplicities. Namely, we will calculate the difference 

A Ql = N Ql - N t , (6) 

1 For generality, we will often use the notation Q in our formulae, having in mind that 
Q means charm quark in the most of the cases (Q = c). 
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analogous to quantity fl3]). Let us note that the multiplicity N L is equal to 
the multiplicity in the e + e~ event taken at the energy W = mw, i.e. Ni = 
NuiW 2 = rriyy), where raw is the W boson mass. 

Hadron multiplicity in a heavy quark event is represented by the following 
equation: 



N Q i(m w , m ) = (n Q + n t ) + N Ql (m w , m ) 



(7) 



with 



N Q i{m 2 w ,m 2 ) = — 

iV 



d 4 q 



m 



w 



x n^(pi,p 2 ,g) 



d A k 



S{{pi+P2 ~ qf ■ 
^ u (q,k,m) n g (k 2 



m 



wj 



(8) 



(2tt) 4 P 

where we have assumed that both bosons are on- shelll Two first terms in 
the r.h.s. of Eq. (171) , are the multiplicities from the leading (anti)quark Q 
and /. They are known from the data. The tensor U^ u (pi,p 2 , q) describes the 
process e + (pi) + e~(j» 2 ) — ► W ± (q) + W T (pi + p 2 — q), where /i, v are the 
Lorentz indices of the W boson with the 4- momentum q (see Fig. [TJ. Two 
blobs in Fig. [Tj are tree diagrams with the j/Z exchange in the s-channel, 
and electronic neutrino exchange in the t-channel. The tensor k,m) 




Figure 1: The diagram describing the process e + (pi) + e~(p 2 ) — > W ± (q) + 
W T (pi +P2~q)i where /i, v are Lorentz indices of the W boson (wavy lines). 
The cut line of the W boson means that it is on-shell particle. 

describes the emission of the gluon jet with the 4-momentum k produced by 
this W boson (see Fig. [2]). The quantity n g (k 2 ) is the average multiplicity 

2 An account of the W boson distribution in its invariant mass q 2 will be discussed 
below (see Eqs. ([50]) - ([53 ]) ). 
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Figure 2: The inclusive distribution of the "massive" gluon jet with the 4- 
momentum k (spiral line). The wavy line is the W boson, whose 4-momentum 
is q, and Lorentz indices are fi, v. 

of the hadrons in the gluon jet with the invariant mass k 2 [TJ [2]. Nq is the 
normalization factor. 

It is convenient to use Lorentz gauge in which the tensor part of the W 
boson propagator has the form: 

dp,{q) = -9nv + ^jf 1 > ( 9 ) 

then 

q^(pi,P2,q) = q u ^( Pl ,P2,q) = . (10) 
The quantity has the following tensor structure: 

$p,(q, k, m) = (-g^q 2 ) C(q 2 , k 2 , qk, m 2 ) + k^k u D(q 2 , k 2 , qk, m 2 ) 

+ (terms ~ q^, q v ) + (term ~ e^ aP q a k P ) . (11) 

Due to condition f TTOT) . terms proportional to q^ or/and g M gives zero contri- 
bution to Nqi after convolution in Lorentz indices with the tensor IL^. 

In the first order in the strong coupling constant, <& ^(q , k , m) is repre- 
sented by the sum of three QCD diagrams presented in Figs. [3j HJ and [5] (the 
crossed diagram is taken with the factor 2). 

One can use the following useful relation 



J i 4 k k^ D(q 2 , k 2 , qk, m 2 ) = (-g^ q 2 ) J d A k 



{qk) 2 - q 2 k 2 



3g 4 

2 7„2 „;„ 



x D(q 2 , k , qk, m ) + (term ~ q^q u ) . (12) 
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Figure 3: The inclusive distribution of the massive gluon jet with the virtu- 
ally k 2 . The wavy line is the W boson, whose 4-momentum is q. The thick 
quark line is a heavy quark, while the thin line is a light quark. The cut 
quark lines mean that these quarks are on-shell quarks. 




Figure 4: The same as in Fig. 10, but with the gluon jet emitted by the light 
quark. 

Detailed analytical calculation^ of the diagrams in Figs. [3]l5] show that 

D = ~{qkf-q^ C ' (13) 

up to small power-like corrections of the type 0(m 2 /m^) and 0(k 2 /m 2 v ). 
The corrections 0(k 2 /m^) can be neglected since it is the region k 2 ~ m 2 
that makes the leading contribution to our main integral (125]) (see Eq. ( 1361) 
and a comment after it). 



3 Explicit analytic expressions for the functions C(q 2 , k 2 , qk, m 2 ) and D(q 2 , k 2 , qk, m 2 ) 
are quite complicated to be shown here. That is why, we will present only the final 



expression ([27] 
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Figure 5: The interference diagram which also contributes to the inclusive 
distribution of the gluon jets with the virtuality k 2 in the W boson. 



The normalization factor N in (jSj) is given by the expression 
1 / m 2 \ 2 / m 2 



.Y : -[} ) ( I + 

7i 



3tt V m 2 w J \ 2m 2 w 



1 4k 6{q2 ~ m2w) 6{{pi +p2 ~ q)2 ~ m2y) { ~ g " Vq2) n ^ 

h J^ 5 (9 2 - m w)S((Pi + P2-q?-ml)(-g^)^. (14) 



Let us define 



Then we derive from ©, and (HTJ-dHfl 

(m w -m) 2 (gfc)_| 



/(iA; 2 C 
-^n g {k 2 ) I d(qk)C(m 2 w ,k 2 ,qk,m 2 ) . (16) 



Qg (ff*0- 
The integration limits for the variable g/c are: 



(gfc)_ = m^v^ ; 

(^) + = (m 2 v + P-m 2 )/2. (17) 



4 Note that the antisymmetric part of $ M „ (i.e. the last term in Eq. (TIT]) ) gives zero 
contribution to Nqi. 
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The hadron multiplicity from the gluon jet with the fixed virtuality k 2 is 
related with the multiplicity from the gluon jet whose virtuality varies up to 

em 



n g {k 2 



k 1 



d 
dk 2 



N 9 (k 2 ) . 



An explicit analytical expression for C(m^, k , qk) shows that 



C{m w , k , qk, m 

qk=(qk) 

After introducing the notation 
G(m^, k 2 , qk, m 2 ) = 
we come to the following equation: 

{raw —m) 2 



C(m w , k , qk, m z 



qk=(qk) A 







N Q i(ra 2 w ,m 2 ) 



dk"_ 



_d_ 



C(q 2 ,k 2 ,qk,m 2 ), 



(18) 



(19) 



(20) 



N g (k 2 ) d(qk)G(m 2 w ,k 2 ,qk,m 2 ) . (21) 



(qk). 



Let us define 



d(qk) G(m w , k ,qk,m ) = Cj 



a s (k 2 



71 



E Ql (m 2 w ,k 2 ,m 2 ) , (22) 



(qk). 



where Cp = (N 2 — 1)/(2N C ), and N c is a number of colors. Eq^iti^, k 2 , m 2 ) 
is the inclusive spectrum of the gluon jet emitted by Q/-quark pair (see 
Ref. [U [2] for more details). Then the multiplicity in heavy quark event (|2"TT) 
can be represented in the form: 



N Q i(m w , m) = (n Q + ni) + C, 



(myy—m) 2 

dk 2 a s (k 2 ) 



N g (k 2 )E Ql (m 2 w ,k 2 ,m 2 " 



k 2 TT 



(23) 

Correspondingly, the hadron multiplicity in the light quark event is given 



by 



NAm 



w) 



2n l + C F j ^^p-N^E^k 2 ), (24) 



where E^m^, k 2 ) = Eq^myy, k 2 , 0). As a result, we obtain the multiplicity 
difference: 

m fdk 2 

A Ql = (ng -n{)- J — N g (k 2 ) AE Ql (k 2 , m 2 ) , (25) 

where the following notation is introduced: 

AE Ql {k 2 , m 2 ) = Et{m 2 w , k 2 ) - E Ql {m 2 w , k 2 , m 2 ) . (26) 

Analytical calculations of the diagrams in Figs. [31 IH and result in the 
following expression: 

G{q 2 , k 2 , qk, 0) - G{q 2 , k 2 , qk, m 2 ) 

a s {k 2 ) 8m 2 qk [k 2 {q 2 - 2qk) + 2m 2 (q 4 - 2q 2 qk + 2{qk) 2 ){qk) 2 } 
F ~ q 2 [k 2 q 2 {q 2 - 2qk) + Am 2 (qk) 2 } ' ( ' 

As before, we have omitted small power-like corrections of the type 0(m 2 / rriyy) 
and 0(k 2 /m 2 w ))E 

One can conclude from Eqs. fl22l . (127]) that AEQi{k 2 ,m 2 ) is a function of 
the dimensional variable 

k 2 

P = — 2 - (28) 

Indeed, let us define: 

qk=[q 2 (l-x) + k 2 -m 2 ]/2. (29) 
The new variable x varies within the limits: 

^ x ^ x max = (l - Vk^/mw^j - m 2 /m 2 w . (30) 
One can safely set x max = 1, and obtain 

. , > 1 /" , . (l-x) 2 (l + x 2 )+Ax 2 p 



= i[2 + p(3p-2)]lni + ~(5 + 6p) 

+ 1 -p(3p-8)J(p) + 6 1 -— M, (31) 
2 p — 4 



5 Remember that we have to put q 2 — in (|27|) . 
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where 



^ln ^5 , p>4, 



J(p) = { 1 , p = 4 , (32) 

/^arctan(^p), p<4. 

Remember that we considered the process e~e + — > W + W~ and com- 
pared two possible subsequent hadronic decays of one of the bosons, W + — > 
cd (cs) + (gluonjets), and W + — > ud (us) + (gluonjets). The quantity 
AE d (k 2 /ml) (|3T|) describes the difference of the distributions of the gluon 
jets in their invariant mass k 2 in these processes (m c is the mass of the charm 
quark). Let us stress that Eq. (J3T1) coincides with the expression for the dif- 
ference of the gluon jet distributions derived in the case when the W + boson 
is a product of the top weak decay t —>■ b + W + (see Eq. (35) from Ref. [5]). 

In terms of variables 

Tfl 2 I 

y = In — = In - , (33) 
£r p 

and 

F = ln^, (34) 
the multiplicity difference looks like 

Y 

A Ql = (n Q -n t )-J dyN g (Y - y) AE Ql (y) , (35) 

— oo 

with Eqi defined by Eqs. f l3Tj) . (1321 . The function AEQi(y) is shown in Fig. [61 
Since 

AE Q (y) -Hg-bl, ( 36 ) 

y— > — oo o 

the integral in (1351) converges rapidly at the lower limit. Asymptotics of 
AEQi(y) at large y is the following: 
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AE Q1 (y) 




Figure 6: The function AEqi(y). 



In Ref. pQ the analogous formula for the multiplicity difference in e + e~ 
events not mediated by the W bosons ([3]) was derived: 



where 



Y 

5 Ql = 2(n Q -m)-J dyN g (Y - y) AE Q (y) 



AE Q {p) = l + pQp-3j ln-+ (l + 7p 



(38) 



+ p (7p - 20) J{p) + 20 , (39) 

with J(p) defined above (1321) . 

The function AEq(y) has the following asymptotics at y — > oo: 

A£ Q (y) ^y-i (40) 

As one can see from Eqs. (HO"]) . (1371) . AEqi(y) = AEq(y)/2 at large y. More- 
over, numerical calculations show that AEqi is very close to AEq/2 at all y 
(see Fig. [7j). Thus, we get the prediction: 
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ae q (y) , 



2AE Q1 (y) 




Figure 7: The function AEq(y) (solid curve) vs. function 2AEqi(u) (dashed 
curve). 

A Ql = ^5 Ql . (41) 
Now we can calculate the total hadron multiplicity from the W boson: 



N w = Ni(m w ) + 







2 + 


Vcs 


2 


\v ud 


2 + 


v us 


2 + \v cd \ 2 + \v cs 


2 



(N d - Ni) 



= Ni{m w ) + -5 d . (42) 

For our numerical estimates, we will use the corrected experimental value of 
5 d from Ref . [1] : 

8 d = 1.03 ± 0.34 (43) 

The light quark multiplicity at the energy W = my/ was recently estimated 
to be 0: 

Niimw) = 19.09 ± 0.18 , (44) 

that results in 

N w = 19.34 ±0.21 . (45) 

Let us estimate effects associated with W boson decays into ub and eft- 
pairs. The account of the 6-quark production increase the multiplicity from 
the W boson (JH by the quantity AN^, i.e. N w -> N w + AA^, where 

AN b w = \V ub \ 2 (N bl - N t ) + \V cb \ 2 (N bc - N t ) . (46) 
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The first term in (14"6"i) is equal to 



2 



V ub \ 2 5 bl . (47) 



In order to estimate the second term, we will use the fact that the emission 
of the massive gluon jets by heavier quark is suppressed, that results in 
E bc < Ebi- Thus, we get the inequality: 

N bc -N t < N bl — Ni — (n c - m) + - 5 bl . (48) 

Using the average values of n c = 2.6, rii = 1.2, and 5 b i = 3.12 [4], as well as 
CKM matrix elements |V u6 | = 3.95 ■ 10" 3 , \V cb \ = 38.6 • 10" 3 [6J we obtain 
from (g6])-(|48D that 

AN^ < 1.2 • 1(T 2 . (49) 

It means that one can ignore the contribution to Nw from the multiplicity 
difference between ub(cb)-eveiat and light quark event. 

Now let us take into account that the W boson with the 4-momenta q in 
Fig. [1] is not on-shell, but has a distribution in its invariant mass q 2 . The 
denominator of the W boson propagator is equal to 

q 2 - mw + im w Tw , (50) 

where Tw is the full W width. As a result, the hadron multiplicity from the 
W boson is given by the formula: 



(W— mw) 



2 



N w -\s d = L j dtf ^—-j-l—pg- F(W\ i) Af,(9 2 ) , (51) 

Ql 



where 



(W—m w ) 



2 



"= I V v-w + ^ W'ft- (52) 



The function F(W 2 , q 2 ) is defined by the invariant part of the diagram in 
Fig. [1] integrated in d 3 q at fixed q 2 . It can be chosen to be dimensionless0 
i.e. dependent on the ratio q 2 /W 2 . 



3 Since g 2 -independent dimensional factors may be safely omitted in F. 
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The main contributions to the integrals in (1511) , (j52p come from the region 
q 2 — mly ~ m w T w . The argument of the function F(q 2 /W 2 ) varies rather 
slowly in this region, contrary to the factor [(q 2 — m^) 2 + m 2 v T 2 v ]~ l which 
has a sharp peak around the point q 2 = m 2 ^. Thus, we can assume that an 
explicit form of F(q 2 /W 2 ) is not important and put F — 1 in Eqs. flSTl) . fl52|) . 
By using a fit of the data on light quark multiplicity Ni(E 2 ), we can estimate a 
possible variation of the multiplicity ( |5T|) . The numerical calculations show 
that at W = 183 GeV (collision energy of LEPII experiments on W + W~ 
production [Z]-[8]) it differs from Ni(m w ) by the quantity 0.050 Then we 
sum the uncertainty 2 x 0.05 = 0.10 with the errors of N w (145]) in quadrature 
and obtain: 

Ni(m w ) = 19.34 ±0.23. (53) 

Our result (153]) is in a good agreement with the experimental data from 
OPAL 0, 

N w = 19.3 ± 0.3 ± 0.3 , (54) 

and DELPHI [8], 

N w = 19.44 ±0.13 ±0.12 . (55) 

Formulae (T4T]) and (14*2]) is our main theoretical result. Let us underline 
that these equations relate measurable quantities, and they do not depend on 
the explicit form of the hadron multiplicity from the gluon jet N g (k 2 ). 

Our main prediction is 

A cl = N cl -N l = 0.52 ±0.17. (56) 

This prediction can be checked by using available LEPII data on hadron 
multiplicities in W + W~ — > qq' lv\ and W + W~ — > qq' qq' events [7], [H] as well 
as future data from the ILC. It will be one more experimental test of the 
universality (i.e. the process-independent character) of the QCD evolution 
in the multiple hadron production. 

Both the multiplicity (153]) . and multiplicity difference (156]) can be also 
measured at the LHC One of the processes to look for mass-dependent effects 
in the QCD evolution is a single W production with its subsequent decay 
into hadrons. Another possibility is a production of two If' s, one of which 
decays in a lepton mode, while the other goes into hadrons with (without) 
charm production. 

7 This difference becomes smaller with the increase of W. 
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3 Conclusions 



In the present paper we have calculated the hadron multiplicity from W 
boson in pQCD. The difference of the hadron multiplicities in the W^-boson 
hadronic decays with and without charm production is also estimated. Previ- 
ously pQ-jS] the multiplicity difference in bb (cc)-event and light quark event 
was calculated on the same ground. The nice agreement of these theoretical 
predictions with the LEP data confirms the universal character of the mech- 
anism of multiple hadron production in QCD via evolution of quark-gluon 
showers. 
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